The hydrino and other unlikely states 
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We discuss the tightly bound (hydrino) solution of theKlein- Gordon equation for the Coulomb potential in 3 dimensions. We 
show that a similar tightly bound state occurs for the Dirac equation for the Coulomb potential in 2 dimensions. These states 
are unphysical since they disappear if the nuclear charge distribution is taken to have an arbitrarily small but non-zero radius. 
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I. THE HYDRINO STATE 

Naudts [1] recently drew attention to the so-called hy- 
drino state of hydrogen, a strongly bound state which has 
a normalisable solution of the Klein-Gordon equation for 
a Coulomb potential in three space dimensions, but no 
analogue in the Schrodinger equation. It has been men- 
tioned by a minority of quantum mechanics textbooks 
in the past [2], which discard it for various ad hoc rea- 
sons: for instance, because it lacks a nonrclativistic ana- 
logue, or because it yields divergent expressions for the 
expectation value of the potential energy. Claims that 
a tightly-bound state of hydrogen exists have, however, 
been widely publicised [3] . Another relativistic equation, 
the Bethe-Salpeter equation, has been shown to have 
anomalous solutions which on closer inspection are un- 
physical [4]. 

In this letter, we show that similar anomalous tightly- 
bound states for a Coulomb potential also exist in the 
Dirac equation in two space dimensions. We claim that 
all these anomalous states are unphysical by developing 
an argument used by Schiff [5], which was first used by 
Loudon [6] for the one dimensional Schrodinger equation 
with the potential —a/ \x\ . 

We shall begin with a brief discussion of the hydrino 
state for the Coulomb potential V(r) = —a/r which is 
taken to be the fourth component of a four-vector. If 
the solution ip(r,9,(p) for the Klein-Gordon equation of 
a particle of mass m is separated out into its radial and 
angular components in the usual way and we consider 
only states of angular momentum I = 0, then if f(r) = rip 
we have (in natural units where H = 1 = c) 
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Consider bound states so that energy E < m. Write 
K = \J m 2 — E 2 , z = 2nr and then 



/"(*) + 



A 

z 



(1/4 -V 2 ) 



m = o, 



(2) 



where A = aE/n and (X = +^/l/4 — a 2 . Eq. (2) is 
called Whittaker's equation [7] and the solutions can be 



expressed in terms of confluent hypergeometric functions. 
To see this take 



/0) = z» +1 / 2 e- z/2 w{z) 



(3) 



to obtain the defining confluent hypergeometric or Kum- 
mer's equation [7] 

zw" + (b~ z)w' - aw = 0; a = fi + 1/2 - A, b = 2fx + 1 

(4) 

The roots t of the indicial equation, giving the lowest 
power of z when w is expressed as a power series, are t = 
and t = l — b= — 2/i. The solutions of Eq.(4) which allow 
the wave function ip to be square integrable to infinity are 
polynomials w{z) — M(a, b, z) with a = —N where N is 
a non-negative integer. 

The conventional choice is t = 0; then Xn = [1 + 1/2 + 
N, whence 



En/tu = 



-1/2 



+ 1/2 + VV4-<* 2 ) 
In particular the ground-state energy Eq is given by 



(5) 



E /m= (1/V2 v/l + y/l - 4a 2 



(6) 



It is easy to show that to lowest order in small a Eq.(5) 
gives the non-relativistic spectrum for the hydrogen atom 
obtained from the Schrodinger equation. 

There are also anomalous solutions of Eq. (2) from 
the other root t — —2[i (or equivalently set fi = 



— \/l/4 — a 2 ). They give the spectrum 
E NA /m= 1 



(A + 1/2- Vi/4-a 2 )' 



-1/2 



and in particular 

E 0A /m -- 
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This is the solution which is called the hydrino. It is 
normalisable but has no counterpart in the Schrodingcr 
equation. Indeed for small a, E a — ma and the electron 
is tightly bound. 

It is tempting to say with Migdal [8] that the one parti- 
cle Klein-Gordon equation cannot be taken seriously and 
field theory is required. But we now demonstrate sim- 
ilar tightly bound states for small a in the one-particle 
Dirac equation, where tight binding should not cause a 
problem [9]. 



II. THE COULOMB POTENTIAL IN THE 
TWO-DIMENSIONAL DIRAC EQUATION 

This problem has been discussed previously [10] and 
anomalous solutions were noticed by Sigurdsson [11]. In 
two space-dimensions we can take the gamma matrices 
7 X , 7 y and 70 to be the Pauli matrices icr X7 ia y and o z 
respectively. Then the Dirac equation for a particle of 
definite energy E in the presence of the potential V(r) = 
—a/r is 



Convert to polar coordinates and take 



ip = 



Mr 9) ( f(r)e i(j - 1/2)e 
~ V ff (r)^' +1 /2)6 



(8) 
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which is an eigenstate of the z-component of total angular 
momentum J z = L z + ^cr z with eigenvalue j. Eq. (8) 
now becomes 



(m - E - a/r) f + dg/dr + (j + 1/2) g/r = 
df/dr - (j - 1/2) f/r +(E + a/r + m)g = 

We follow the approach of Tutik [12]. Write 

f = Fy/(l + E/m), g = Gy/{\ - E/m) 

and $^2 = F ± G, and scale to z to obtain 

d$i $1 .$2 , U a 

— — + - j h -$1 = — {E$! + TO$ 2 ) 

dz 2z z 2 nz 
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Take $1,2 = ^1,2 exp(— z/2); then 



(11a) 
(lib) 

(12) 



so (f>2 can be written in terms of </>i. The second order 
differential equation for is 
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or 



Write 



zcp'{(z) + (2 - z)^ - (1/2 - oE/k) 4>i 
(f-l/4-a 2 )fa/z 

1 = z'w so we obtain 



zw"(z) + {2 + 2t- z)w' + (t 2 + t)w/z - tw 
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Hence take 



or 



to obtain 



t 2 +t 



1/4 - a 2 



t = -1/2 ± y/j 2 -a 2 



(17) 
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which again is Kummer's equation. The regular solution 
M(a, b, z) is square integrablc to infinity only if it reduces 
to a polynomial, which it does if 



a = t + 1/2- aE/K = —N 



(19) 



where N is a non-negative integer. 

Conventionally one chooses the positive sign in Eq. 
(17) which yields 



aE N / K = N + v/j 2 - a 2 



or for j = \ 



En 
m 



1 + 



a 



(N + y/l/A -a 2 ) 2 
and the ground state is given by 



-1/2 



E Q /m = a/1 - 4a 2 



(20) 



(21) 



(22) 



As r — > 0, /, g ~ r* = r5 + Vr^ so that the wave func- 
tions remain square-integrable down to the origin even 
though they diverge. They are the relativistic generali- 
sation of the solutions of the Schrodingcr equation in two 
dimensions which has the energy spectrum given by Eq. 
(21) to lowest order in a 2 . 
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The anomalous solutions arise from the other root of 
Eq. (17): t = -1/2 - ^1/4- a 2 . They too are square 
integrable down to the origin. The anomalous energy 
eigenvalues are given bya = t+ l/2 — aE^/n = — N or 



aE NA /n = N- Vl/4- a 2 



(23) 



Because k > 0, the lowest anomalous state with N = 
satisfies 



whence 



aE 0A /n = a 2 



EoA/m = -y/l - Aa 2 



(24) 



This solution has no non-relativistic analogue and like the 
hydrino corresponds to a state which is tightly bound for 
small a. 



III. DISCUSSION 

We believe that the three dimensional Klein-Gordon 
hydrino state and its two dimensional Dirac analogue 
are unphysical in spite of having normalisable wave func- 
tions. Their obvious failings are that 

1. they lack non-relativistic counterparts even for ar- 
bitrarily small coupling; 

2. the states persist even when the coupling is turned 
off; 

3. the strength of the binding increases as the coupling 
strength a decreases. The maximum binding oc- 
curs for a — when the potential has disappeared 
completely 1 . 

These reasons may be sound but they are not deci- 
sive either mathematically or physically. We now demon- 
strate that, if the point charge of the nucleus is replaced 
by a charge extending over an arbitrarily small but finite 
radius R, then the anomalous functions become unac- 
ceptable because for small enough R they cease to satisfy 
the appropriate wave equation. We will prove this for the 
Klein-Gordon equation; the argument for the Dirac case 
runs similarly. 

For simplicity we look at a conducting spherical shell 
of radius r = R, so that V(r < R) = —a/R while V(r > 
R) = —a/r. We need consider only a 2 < 1 since we are 
investigating the bound-state energy in the limit a — > 0. 
Then the inner solution of Eq. (1) is 

f(r < R) = Csm(qr), q = ^{E + a/R) 2 - m 2 (25) 



1 We could call these anomalous states homeopathic states 
because the smaller the coupling, the larger the effect. 



For R small enough to make a/R » m, one has q — > 
a/R, so that qr = ar/R < a < 1 f(r) = Cqr(l + 
0{{qr) 2 ) . In terms of the scaled variable z, and defining 
Z = 2nR, this entails 



zdf(z)/dz 



1 + O(Z) 



(26) 



This solution has to be joined to the outer solution of 
Eq. (2) for z > Z = 2kR. For w(z) defined by (??) the 
general solution of Eq. (4) reads [7] 

Bz^Mia- b+ 1,2 -b,z) (27) 



w = AM(a,b,z 
where asz^oo 



M(c, d, z — > oo ) 



T(d) 



exp(z)z 
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Hence 

w(z - 
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A 



r(c) 

r(6) r(2 - b) 

T(a) r(a-6+l) 



(28) 



exp(z)^ 



a—b 



(29) 

For a bound state the expression in the brackets must 
vanish so that 



w = M(a, b, z) — 

T{b)T{a-b+l) j 



r(a)r(2-6) 
Thus for small z 

w(z -> 0) -» 1 - 

r(6)r(a-6 + 



z 1 ~ b M(a-b+ 1,2 -b,z) (30) 



1) 



r(a)r(2 - b) 



.1-6 



(a-b+1) 
(2-6) ; 



2-b 



(31) 



hence 



/(«->())-> (1- z/2 + ...)z 



p+l/2 



W 



(32) 



Recall a = jj, + 1/2 - aE/n, b = 2fi + 1, \i = - a 2 , 



2-b ~ 2a 1 . 



so for small a 

l-6~ -l + 2a 2 , 

Since we are considering the limit as a — > 0, we can 
disregard a 2 in the various exponents of z and so as z — > 



/ 



r(6)r(o- b + 1) 



(l + z + ...) 



(33) 



r(o)r(2-6) 

Hence requiring f(z) to be proportional to z for small z 
in order to satisfy Eq. (26)gives 

l/r(a) = => a = -N, AT = 0,1, 2,... (34) 

This is precisely the conventional quantisation condition 
which gives the energy eigenvalue spectrum of Eq. (5) 
thus eliminating the anomalous solutions. 

For the two dimensional Dirac equation for a Coulomb 
potential a similar argument, with Bcssel functions re- 
placing trigonometric functions for the inner solution, 
will give f(r = R) — > const; g(r = R) — > for R small 
enough and hence <pi — z t w — > const. This is only pos- 
sible if the positive sign is taken in Eq. (17), ruling out 
the anomalous solutions here too. 
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IV. CONCLUSIONS 

We have shown that there are normalisable solutions of 
the Dirac equation for the Coulomb potential in two di- 
mensions which are tightly bound, just like the hydrino 
solution of the Klein-Gordon equation in three dimen- 
sions. We have displayed the unphysical nature of such 
solutions by smearing the nuclear charge distribution so 
that it has an arbitrarily small non-zero radius; in that 
case there are no anomalous solutions. We suggest that 
outside of science fiction this is sufficient reason to dis- 
regard them. We leave to others to find a general crite- 
rion to exclude the anomalous solutions for the Coulomb 
potential with a point charge in two and three space di- 
mensions: in the Appendix we indicate a possible route 
to this end. 
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V. APPENDIX 

The requirement that wave functions have to be square 
integrable is not the only restriction on physical states in 
quantum mechanics. As de Castro has pointed out in his 
study of Klein Gordon particles [13] for a similar prob- 
lem, it is also necessary to ensure that the Hamiltonian 
is Hcrmitian. We shall not attempt to do this here, but 
a necessary consequence of Hermiticity will be that the 
eigenfunctions of the Hamiltonian corresponding to two 
distinct energy eigenvalues are orthogonal. 

The eigenfunctions fo(z) for the normal ground state 
energy given by Eq. (6) and Joa(z) for the hydrino energy 
given by Eq. (7) are similar. They are both of the form 
given by Eq. (3): 

f(z) = z^ +1/2 e~ z/2 w(z) (35) 

where in both cases w(z) = M(a, b, z) and a = 0; 
i.e. w(z) is constant. Neither fo{z) nor /oa{z) van- 
ishes except at z = 0. Hence the hydrino wave function 
ipoA = foA(z)/r cannot be orthogonal to tpo = fo(z)/r. 
Thus the hydrino state cannot be physical provided that 
it is possible to ensure that the Hamiltonian for the hy- 
drogen atom in the Klein-Gordon equation is Hcrmitian. 
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